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ABSTRACT 

We demonstrate that tidal evolution of the inner planet (‘e’) of the system orbiting 
the star p 55 Cancri could have led to passage through two secular resonances with 
other planets in the system. The consequence of this evolution is excitation of both 
the planetary eccentricity and inclination relative to the original orbital plane. The 
large mass ratio between the innermost planet and the others means that these exci¬ 
tations can be of substantial amplitude and can have dramatic consequences for the 
system organisation. Such evolution can potentially explain the large observed mutual 
inclination between the innermost and outermost planets in the system, and implies 
that tidal heating could have substantially modified the structure of planet e, and pos¬ 
sibly reduced its mass by Roche lobe overflow. Similar inner secular resonances may 
be found in many multiple planet systems and suggest that many of the innermost 
planets in these systems could have suffered similar evolutions. 

Key words: planet-star interactions - planets and satellites: dynamical evolution 
and stability 


1 INTRODUCTION 

The importance of secular gravitational interactions in 
sculpting planetary systems has long been appreciated as 
being necessary to understand the regularities of our own 
solar system. It is known that the long term stability of the 
terrestrial planets of our solar system is imperilled by the 
influence of secular gravitational interactions between Mer¬ 
cury, Venus and Jupiter (Laskar 1994; Batygin & Laughlin 
2008; Laskar & Gastineau 2009). With the discovery of a 
host of extrasolar planet systems over the past twenty years, 
including many multiple systems (Wright et al. 2009; Lis- 
sauer et al. 2011a; Fabrycky et al. 2014), the role of secular 
interactions in determining the stability and architecture of 
planetary systems is becoming of ever greater importance. 

The planetary system around the star p 55 Cancri is of 
particular importance in this regard. One of the first sys¬ 
tems discovered by radial velocities (Butler et al. 1997), it 
also provided one of the first extrasolar giant planets on or¬ 
bital scales similar to that of Jupiter (Marcy et al. 2002) and, 
most recently, the revision of the innermost planet parame¬ 
ters (Dawson & Fabrycky 2010) resulted in the photometric 
detection of a planetary transit (Winn et al. 2011; Demory 
et al. 2011) which provides some of the best constraints on 
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the properties of rocky super-earths (Gillon et al. 2012; De¬ 
mory et al. 2012). Furthermore, the presence of a Jovian 
class planet with an orbital period of 14 days (55 Cane b) 
places this system in the enigmatic class of ‘warm Jupiters’, 
which don’t fit naturally into the evolutionary paradigms 
developed to explain the presence of Jovian planets with 
shorter periods. Finally, reports of a tentative detection of 
an astrometric signal with the HST Fine Guidance sensor 
(McArthur et al. 2004) suggest that this planetary system 
may possess a substantial dispersion in orbital inclinations, 
with a nominal offset of 37° ± 7° between the outer planet 
55 Cane d and the edge-on inner planet 55 Cane e. 

It has also been suggested (Kaib, Raymond & Duncan 
2011) that the presence of a distant M-dwarf companion 
(Mugrauer et al. 2006) can induce a coherent precession of 
the entire planetary system (Innanen et al. 1997), potentially 
leading to a substantial obliquity offset between the stellar 
spin and the orbital plane of the planetary system. However, 
for the 55 Cancri system, spin-orbit coupling between the 
star and the planets likely quenches such an offset unless the 
companion is in a highly eccentric (e ~ 0.95) orbit (Boue’ & 
Fabrycky 2014). Observationally, the situation is unclear as 
Bourier & Henrard (2014) claim to have detected a projected 
obliquity of 72±12°, which Lopez-Morales et al. (2014) claim 
is impossible given the very low stellar rotation. 

All of this recent activity suggests that the 55 Can- 
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cri system is interesting for both its dynamical architecture 
and history as well as for the composition of the component 
planets, in particular 55 Cane e, as an example of a mas¬ 
sive, earth-like planet transitting a bright star. Our analysis 
will focus on the fact that 55 Cane e is close enough to the 
star that it should have experienced substantial tidal dis¬ 
sipation, and consequent orbital evolution. Such evolution 
can have an effect on more distant planets too, as secular 
coupling between the planets can transfer angular momen¬ 
tum back and forth between the planets in the system (e.g. 
van Laerhoven & Greenberg 2012). In § [2] we will discuss 
our model for tidal interactions and describe the evolution 
of the system within the context of classical secular theory, 
and we will show how this can provide a potential explana¬ 
tion for the large inclination dispersion. In §[3] we investigate 
the evolution further using numerical integrations, to clarify 
the dynamical limitations on the excitation of inclination 
and discuss the influence of the near-resonance of Planets 
55 Cane b and 55 Cane c on the secular oscillations of the 
system. Finally, we discuss the possible influence of tidal 
heating on the structure of 55 Cane e and on the likelihood 
of other known planetary systems where inclination excita¬ 
tion could have occurred. 


2 TIDAL EVOLUTION IN THE CLASSICAL 

SECULAR APPROXIMATION 

The physical parameters of the innermost planet, 55 Can- 
cri e, are estimated to be a mass of 7.99 ± O.25M0 and a 
radius of 2.OO±O.14R0 (Nelson et al. 2014; Winn et al. 2011; 
Demory et al. 2011). The radius of a perovskite planet of this 
mass is ~ l.Oi?®, using the fitting formulae of Seager et al. 
(2007). Thus, although it is closer in mass to a Neptune or 
Uranus, this is expected to be a predominantly rocky planet 
(although it could have a finite water contribution as well). 
We wish to therefore consider a range of tidal dissipation 
extending from terrestrial values to ice giant values. 

To implement the tidal dissipation within the secular 
formalism, we use the same approach as in Hansen & Mur¬ 
ray (2014), following on from prior work of Wu & Goldreich 
(2002) and Greenberg & van Laerhoven (2011). We have 
also previously examined tidal dissipation for extrasolar gi¬ 
ant gas planets (Hansen 2010), using the formalism of Hut 
(1981) and Eggleton et al. (1998). However, the restriction 
of classical secular theory to the quadratic level in eccentric¬ 
ity means that we will restrict our tidal expressions to the 
same order. We must also choose a dissipation level more ap¬ 
propriate to terrestrial planets. Choosing a bulk dissipation 
constant that yields Q' = 10 for the Earth (e.g. Goldreich 
& Soter 1966), we find 
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We have included the estimated mass and radius for 
55 Cane e in these estimates, but let the semi-major axis 
float. This demonstrates that we expect substantial tidal 
evolution for a planet like this even out beyond 0.1 AU, 


unless the strength of the dissipation is several orders of 
magnitude weaker than expected. 

Of course, there is little evolution if the orbit starts 
close to circular, but the presence of much larger planets in 
close proximity suggests that a finite eccentricity would be 
excited by secular perturbations anyway, resulting in tidal 
evolution (e.g. Mardling & Lin 2004). Thus, we need to con¬ 
sider the fact that the currently observed configuration could 
be the result of substantial tidal migration if we wish to infer 
something about the origin and evolution of this planetary 
system. 

2.1 Secular Architecture 

Consider then the classical secular solution for the 55 Can- 
cri system with the currently measured parameters of Nel¬ 
son et al. (2014). In addition to the classical equations to 
second order in eccentricity and inclination (Murray & Der- 
mott 1999), we also include the effects of relativistic pre¬ 
cession (Adams & Laughlin 2006). There is potentially also 
an additional precession due to the presence of a distant M 
dwarf companion, a 0.26Mq M dwarf located at 1260 AU 
(Mugrauer et al. 2006). However, this is too small to affect 
the dynamics of the planetary conhguration on small scales, 
although it can generate precession of the system as a unit 
(Innanen et al. 1997; Kaib et al. 2011; Bone’ & Fabrycky 
2014). In order to better match the observations, we fur¬ 
thermore assume that the planetary system invariable plane 
is inclined at 37° to the line of sight, except for the transit- 
ting 55 Cane e. This increases all the masses by a factor of 
1.25 from that inferred in the radial velocity solution. These 
are summarised in Table [T] 

The modal structure for the eigenvalues that results 
from this configuration is such that the eccentricities of the 
outer two planets are decoupled from the inner three (as 
noted by van Laerhoven & Greenberg 2012), whose interac¬ 
tion will dominate the eccentricity of 55 Cane e in the ab¬ 
sence of dissipation. On the other hand, all four of the inner 
planets are coupled to a low frequency inclination eigenmode 
that represents the precession of the inner system under the 
influence of the M dwarf perturber and 55 Cane d. 

If the inner planet migrated due to tidal dissipation, 
then the secular structure of the system will have changed 
as well. Figure [T] shows how the secular eigenfrequencies of 
this system vary as we vary the semi-major axis of Planet e, 
keeping all masses and other semi-major axes fixed, from 
the current location to a semi-major axis of 0.05 AU, cor¬ 
responding to a period of 4.2 days. The two lowest eigen¬ 
frequencies in eccentricity correspond to the outer planets d 
and f, and do not change because the outer planets are es¬ 
sentially decoupled. The lowest frequency inclination eigen¬ 
function corresponds to a mode coupling the four innermost 
planets. 

The higher frequency eigenmodes show a more inter¬ 
esting structure. In particular, we see that the two high¬ 
est frequency modes in both inclination and eccentricity 
experience resonant behaviour in the interval Oe ~ 0.026- 
0.031 AU. Thus, if 55 Cane e starts with a semi-major axis 
exterior to 0.032 AU (orbital period of 2.15 days), it will 
cross first an inclination resonance and, shortly thereafter, 
an eccentricity resonance. These will result in mode mixing 
and can result in interesting behaviour, as we shall see. For 


© 0000 RAS, MNRAS 000, 000-000 












Tilting Orbits 3 


Table 1. Basic Architecture of the 55 Cancri Planetary System, based on the parameters of Nelson et al. 
(2014) and assuming an inclination of 37° to the line of sight for all but the innermost planet. 


Planet 

a (AU) 

M (Mj) 

e 

0.035 

0.0251 

b 

0.1134 

1.056 

c 

0.2374 

0.223 

f 

0.774 

0.184 

d 

5.451 

4.79 



Figure 1. The black curves show how the precession frequencies 
of the secular eccentricity eigenvalues gi change as the semi-major 
axis of 55 Cane e is varied, with all other semi-major axes and 
masses held constant. The green curves show the equivalent ab¬ 
solute values of the frequencies for the inclination eigenvalues, 
|/|. The inset shows a zoom around the most interesting part of 
parameter space. The various curves are labelled with the sym¬ 
bols for all planets that make at least a 10% contribution to the 
unit eigenvector of that mode. We see that, if 55 Cane e moves 
through the region ~ 0.026—0.031 AU, the system will experi¬ 
ence secular resonances between both eccentricity and inclination 
eigenmodes. In particular, these highlight that the principal inter¬ 
action is between 55 Cane e and 55 Cane c, not the second most 
distant 55 Cane b. However, all three planets interact via a third 
eccentricity mode, while the outer two planets are essentially de¬ 
coupled. In the case of the inclinations, the inner four planets also 
precess as a unit under the influence of the outermost planet. For 
completeness, the dashed curve shows the value of the relativistic 
precession frequency at the location of 55 Cane e. 

the purposes of the following discussion, we designate the 
eccentricity resonance as gee and the inclination resonance 
as /ce- 

The dashed line in Figure [T] also shows the strength 
of relativistic precession on these close planetary scales. We 
have repeated the above calculation without either of the ex¬ 
ternally imposed precessions and verified that the behaviour 
is qualitatively similar. The relativistic contribution has two 
quantitative consequences. The first is that the frequencies 
of the inclination and eccentricity precessions for 55 Cane e 


are the same on small scales without relativity (the black- 
green pair labelled e in the inset lie on top of one another 
for fle < 0.026AU) and the turnup in the eccentricity eigen- 
frequency of 55 Cane e at Ue < 0.018AU is a consequence of 
relativity. 


2.2 Tidal Evolution 

Let US now consider the tidal evolution of a 55 Cancri ana¬ 
logue in which 55 Cane e starts with an original semi-major 
axis ae(0) > 0.04 AU. The secular coupling of the inner three 
planets means that the tidal dissipation will remove energy 
from the system as a whole instead of just 55 Cane e (Wu & 
Goldreich 2002; Greenberg & van Laerhoven 2011; Hansen & 
Murray 2014), and so, the current rather small eccentricities 
of 55 Cane b and c may have been substantially damped. 
We take the eccentricity of 55 Cane f to represent a pri¬ 
mordial value because it is largely decoupled from the inner 
three planets. Let us start then, for illustrative purposes, 
with the same eccentricity for planets b, c and f, ~ 0.32. 
Figure [2] shows the evolution of the system under the influ¬ 
ence of tidal dissipation using equations ([T]) and ([2]). The 
modelling of the interactions between the tides and the sec¬ 
ular oscillations is performed using the formalism of Hansen 
& Murray (2014). Figure [2] shows three sets of curves, for 
three different levels of dissipation Qg = 10, 100 and 1000. 
We see that the consequence of the resonant crossings is the 
excitation of inclination and eccentricity for 55 Cane e. The 
excitation of the eccentricity also results in a rapid increase 
in the strength of the tidal dissipation and a sharp decrease 
in semi-major axis. 

The size of the rapid ’jump’ in semi-major axis is de¬ 
pendant on the strength of dissipation Qg, and is larger for 
weaker dissipation. This is a little counter-intuitive but can 
be understood in terms of the speed at which 55 Cane e 
moves through the resonance. The origin of the high ec¬ 
centricity of 55 Cane e is the mixing of the two mode ampli¬ 
tudes during the Qce passage, and a slower passage results in 
a larger transfer of amplitude between the two modes. The 
higher eccentricity means a smaller periastron and more in¬ 
ward migration. 

The increase in inclination fg seen in Figure [2] offers an 
explanation for the enrions misalignments of 55 Cancri e and 
d. If 55 Cane e began in the same disk plane as d and all 
the other planets, but migrated inwards through the secu¬ 
lar resonances, then inclination offsets ~ 20° or greater are 
quite plausible. 
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Figure 2. The upper panel shows the temporal evolution of the 
semi-major axis of 55 Cane e as the planetary system evolves 
under the influence of tides and secular coupling. We show the 
evolution for three values of Qg=10, 100 and 1000. We see that 
weaker dissipation actually leads to more orbital migration. The 
reason for this counter-intuitive behaviour can be found in the 
middle panel, which shows that the eccentricity excitation during 
the resonant crossing is stronger if the tidal dissipation is weaker. 
The lower panel shows the increase in the planetary inclination 
due to the crossing of the inclination resonance. 


2.3 Analytic Estimates 

Minton & Malhotra (2011) have investigated the excitation 
of eccentricity in solar system asteroids by sweeping secular 
resonances driven by outer planet migration. This formalism 
can be adapted to understand the excitation of eccentricity 
and inclination in the 55 Cancri system too. One important 
difference is that the interaction discussed in MM 11 is driven 
by the change in the frequency of the external resonance, 
whereas the interaction here is driven by the inward drift 
of the excited body, as it crosses a fixed external resonance. 
In appendix 1X1 we discuss the MMll analysis, demonstrate 
that the end result is similar despite the modification to 
the model, and present the equivalent expression for the 
inclination. Here we will briefly review the outline of those 
results. 

MMll start with a Hamiltonian of the form 

He = —goJ + e'/Zj cos(tDp — oj) (3) 

where J = \/a{i — %/l — e^), go is the precession frequency 
for the longitude of periastron oj, and e parameterises the 
strength of the coupling to the externally maintained secular 
mode, whose orientation is ojp. To this Hamiltonian, we add 
the terms that describe the corresponding interaction for the 
inclination 

Hi = —hoZ + 5\/2Z cos(flp — Q.) (4) 

and Z = a{l — e^)(l — cosi). 



Figure 3. The dashed lines enclose the range of possible final 
eccentricities, depending on initial eccentricity and longitude of 
periastron, calculated in the limit of low eccentricity. The points 
show the solution using the full expression for J, which preserves 
the basic structure of the low e expansion, but shifts the solution 
to higher values. These solutions are for the case = 0.1, and 
will shift according to this quantity and the strength of the tidal 
dissipation. 


For the size of the eccentricity jump from Ji to J/, they 
derive 


J f — Ji -\- 


2|s-o| 


+ e. 


27TJi 

-r-r-r cos OJi. 
l5o| 


(5) 


The size of the jump is thus regulated by the rate of change 
of go due to tidal dissipation, as seen in the previous section, 
as well as the phase of the pericenter precession, Cji. The 
expression for this is 


n a I rrij 2 15 GMc 

^ 8 a I ^ ~ T c^a 
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where n is the orbital frequency, Oj = a/fflj is the ratio 
of semi-major axes of 55 Cane e and perturber j, 7 ^ = 


^^72 + 


- 2ab^}l + b?'> 


and the functions b^s'* 


5/2 ^ '=5/2 

the Lagrange functions (functions of aj). 

Taking a/a from equation o we see that the rate of dis¬ 
sipation actually depends on the eccentricity, which is the 
quantity being excited during the transition. Thus, the speed 
at which the planet crosses the resonance accelerates as the 
eccentricity is excited - a naturally self-limiting process. To 
estimate the final eccentricity, we thus express go = e^go/ro, 
where e/ is the final eccentricity. In the limit of low eccen¬ 
tricity, J ~ ^y/ae^, so that the final expression is 
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e/ 
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Ci + 


we^To 
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2'KTo Bi cos cu 
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(7) 


The size of the eccentricity jump is large enough, how¬ 
ever, that we should use the full expression for J/. Figure[3] 
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shows the resulting range of final eccentricities as a function 
of initial eccentricity, for our nominal Q'^ ~ 10, assuming 
Cb = 0.1. We can see that, for small initial eccentricities, 
there is a unique final eccentricity ~ 0.56, although this 
will, of course, be damped further by tidal evolution after 
the resonance has been crossed. In the absence of further 
eccentricity excitation, this also predicts a final semi-major 
axis ~ ao(l — e^) ~ 0.0205AU, for ao = 0.03 AU. The 
size of the eccentricity jump also scales with the eccentric¬ 
ity of 55 Cane b. If Cb is as small as 0.01, then the final 
eccentricity is ~ 0.18 in the limit of low initial eccentricity. 
Nevertheless, we see that substantial eccentricities can be 
generated even if the eccentricity of the perturbing planets 
are relatively small. The eccentricity of Planet b may also 
experience moderate pumping due to the coupling between 
planets e, b and c via a third mode (see Figure [T]), but this 
is limited by the large mass ratio. 

We find a similar inclination excitation, 

Zf = Zi + OI^+S O^cos^i. (8) 

2 h-o \/ /lo 
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T ^ T 
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Since to once again depends on the eccentricity, we note that 
the inclination resonance lies slightly outside the eccentricity 
resonance, so that we expect this to scale with a, not e/. In 
the limit of low eccentricity and inclination, the analogue to 
equation 0 is 


•2 , 

— + 


2/io%/aef 




27rT'o n 
ho^/a Ci 


cos fli. 


( 9 ) 


In the limit of low initial eccentricity and inclination, this 
still predicts a final inclination ~ 45°, for values of 55 Cane b 
of Ci, = 0.1 and ib = 2°. 

Therefore, this analysis supports the idea, expressed in 
the previous section, that substantial changes in inclination 
and eccentricity are possible for reasonable planetary pa¬ 
rameters. It also demonstrates the dependance on param¬ 
eters, in that eccentricity and inclination changes increase 
with the dissipation timescale tq, i.e. crossing the resonance 
more slowly increases the excitation. 


2.4 Additional Planets 

The present analysis is based on the secular behaviour of 
the observed planets. However, the secular oscillations de¬ 
pend on the effects of all the planets in the system, which 
raises the question whether it is possible to substantially 
alter the properties of the system evolution if there are as- 
yet-undetected planets present. 

To assess this effect, we have taken the collection of 
extant radial velocity data from Nelson et al. (2014), Endl 
et al. (2012) and Fischer et al. (2008) and determined how 
large a planet could remain undiscovered in this data. The 
orbital periods of the five known planets were held fixed at 
the values determined by Nelson et al. but other parameters 
such as mass, eccentricity and longitude of periastron were 
allowed to vary, and then the data was refit with a series of 
models with an additional planet with a range of mass and 
semi-major axis. Orbital periods from 94 to 3900 days were 
searched, and the mass constraints were defined as those 
that would yield a signature larger than a 90% confidence 


Figure 4. The solid points show the currently known planets 
of the p 55 Cancri system, with the radial velocity values for the 
masses shown. The upper limits indicate the 90% confidence mass 
limits for selected regions of period space. Planets with this pe¬ 
riod and masses larger than this limit should have already been 
detected in the radial velocity data. The cross shows the approx¬ 
imate location of the tentative 60.76 day signal discussed in pass¬ 
ing by Nelson et al. (2014). 


level fluctuation with respect to the best fit. The resulting 
mass limits are shown in Figure 01 

For planets in the range 0.4-0.75 AU, radial veloc¬ 
ity signals would be detected for planets with masses > 
4 X 10“®Mq. If we include a single planet of this mass and 
period in our secular evolution, it results in no qualitative 
changes to the evolution. This is because the extra planet 
only exhibits a substantial coupling to two modes. The first 
is a mode dominated by the new planet alone, and the second 
mode is the mode coupling all of the inner planets together 
(the mode marked e-b-c in Figure [TJ. Thus, another planet 
in this gap would share in the mode that mixes all the plan¬ 
ets together, and can contribute to the angular momentum 
budget of the secularly coupled system, but it does not sub¬ 
stantially alter the structure of the modes that undergo the 
Qce and /ce resonant interactions). We have tested the ef¬ 
fects of planets an order of magnitude larger than our RV 
limit and still find no substantial qualitative change in the 
mode behaviour. Nelson et al (2014) note the presence of a 
weak signal in their data with period 60.76 days and velocity 
semi-amplitude ~ 2m/s. This is of similar magnitude to the 
limits we find, and so, if such a planet exists, it would not 
substantially change the secular structure of the system. 

We have also tested the effects of a planet in the gap 
between planets f and d. If such a planet improves the sec¬ 
ular coupling of the inner planets and the outer planets, 
it might change the evolution of the system by providing 
a greater reservoir to store angular momentum and might 
help to drive an even stronger evolution. However, with mass 
limits < O.lMj, extra planets with a > 1 AU appear to only 
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weaken any residual coupling between inner and outer plan¬ 
ets. New planets interior to this would be within the 3:2 
commensurability with planet f and should provide strong 
dynamical signatures that would be detected in the radial 
velocities. 

We therefore conclude that the behaviour discussed here 
is likely to be robust with respect to the presence of as-yet- 
undetected planets in the 55 Cancri system. 


3 NUMERICAL EXPLORATION 

The analytic and semi-analytic results offer a qualitative ex¬ 
planation for the possible origin of a large inclination differ¬ 
ence in the 55 Cancri system, but there are several reasons to 
be cautious about the quantitative results. Our initial esti¬ 
mates suggest that eccentricities and inclinations can reach 
substantial values, exceeding the approximations of classi¬ 
cal secular theory. Furthermore, the proximity of planets 
55 Cane b and 55 Cane c to the 3:1 mean motion resonance 
could possibly amend the secular oscillations and shift the 
resonance locations. This can be incorporated into the clas¬ 
sical formalism in the case of first order resonances (Mal- 
hotra et al. 1989) but higher order resonances scale with 
higher powers of eccentricity and inclination, which would 
introduce nonlinearities into the classical secular equations. 
Thus, we have performed numerical integrations of model 
systems with the MercuryG code (Chambers 1999) to exam¬ 
ine the secular structure of actual 55 Cancri analogues. 

The standard MercuryG algorithm does not account for 
relativistic precession nor tidal evolution of the orbit. We 
have thus included a central potential contribution to pro¬ 
duce the former and account for the latter by instituting a 
drag force by including a radial velocity damping term. The 
small stepsizes required to resolve these short period orbits 
means that we are able to reproduce the expected preces¬ 
sion rate without needing to incorporate invididual stepsizes 
(e.g. Saha & Tremaine 1994). We describe the calibration of 
the latter in appendix [B] 


3.1 Test Particle Orbits 

The simplest test of the effect of the 3:1 resonance is to 
examine the secular oscillations of test particles interior to 
55 Cane b. We simulated several analogues of the 55 Cane 
planetary system with 55 Cane e replaced by a set of test 
particles on co-planar circular orbits and with a range of 
semi-major axes, in order to understand the secular forcing 
this planet would experience from the other four. We chose 
the longitudes of periastron and ascending node randomly 
for the four planets and integrated several cases, to under¬ 
stand what sort of secular behaviours might be expected. 

FigureO shows the expected result, based on the consid¬ 
erations of § [2l We start with 800 test particles on circular, 
coplanar orbits between 0.015 and 0.065 AU, perturbed by 
the four planets of the 55 Cancri system with a> 0.1 AU. 
The figure shows the resulting forced eccentricity and incli¬ 
nation after 1.1 Myr of evolution. Over most of this range, 
the appearance of simple Lagrange oscillations are evident, 
but at a ~ 0.03AU, both eccentricity and inclination grow 
due to the resonance of the free Lagrangian precession of 



0.02 0.03 0.04 0.05 0.06 

a (AU) 


Figure 5. The upper panel shows the eccentricity, while the 
lower panel shows the inclination, after 1.1 Myr of evolution. Fur¬ 
ther evolution can drive the eccentricity to high enough values to 
impact the star. 

the test particles with the secular oscillations of the forcing 
planets. 

However, other behaviour is also possible. Figure [6] 
shows the evolution of the same test particle setup but forced 
by a planetary system with a different choice of eccentric¬ 
ities and longitudes of periastron and ascending node. In 
this case the inclination resonance occurs at approximately 
the same location, but the eccentricity resonance occurs at 
a ~ 0.05AU. Further examination of the results from this 
integration shows that the resonant angles generated by the 
pair 55 Cane b and 55 Cane c circulate most of the time, 
but undergo intermittent passages through resonance, which 
occur when the secular oscillations bring one of the plane¬ 
tary eccentricities close to zero. The origin of this behaviour 
is discussed in more detail in appendix o but the salient 
feature for this analysis is that this interaction between res¬ 
onant and secular behaviour is enough to shift the secular 
resonance in some cases. 

3.2 Tidal Evolution 

In order to examine the consequences of this kind of evolu¬ 
tion for plausible 55 Cancri origins, we perform the following 
numerical experiment on the 55 Cancri system, using the 
MercuryG code. 

We integrated 100 systems with 55 Cane e originally 
located at 0.033 AU (i.e. just outside the secular resonance 
location in the classical estimate). The other four planets 
assume the semi-major axes determined by Nelson et al. 
(2014) from the latest compilation of the data. We assume 
masses for the four outer planets assuming an inclination 
angle of 53° (i.e. an enhancement by a factor of 1.25 above 
the nominal edge-on estimate), although we retain the nom¬ 
inal mass of 55 Cane e because it is transiting and therefore 
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7:2 



0.02 0.03 0.04 0.05 0.06 

a (AU) 


Figure 6. The upper panel shows the eccentricity, while the lower 
panel shows the inclination, after 0.3 Myr of evolution. In this case 
we see that perturbations due to the proximity of the 3:1 reso¬ 
nance have shifted the location of the eccentricity resonance with¬ 
out shifting the location of the inclination resonance markedly. 
The larger eccentricities also lead to enhanced contributions from 
mean motion resonances with 55 Cane b as well. 


has inclination of 90°. Tidal dissipation tends to align the 
longitudes of pericenter, so we draw the initial longitudes of 
perihelion and ascending node, as well as the mean longitude 
at the start, from uniform distributions. The inclinations are 
chosen from a Gaussian distribution with dispersion of 2°, 
which is consistent with estimates for the flatness of com¬ 
pact planetary systems based on Kepler studies (Lissauer et 
al. 2011a; Tremaine & Dong 2012; Fang & Margot 2012). 
Tremaine & Dong (2012) note that larger inclination dis¬ 
persions are possible, but our assumption represents a con¬ 
servative starting point in order to illuminate the size of 
the effect produced by the secular resonance. Although the 
55 Cancri planets have measured eccentricities, these could 
have been substantially reduced by the effects of tidal damp¬ 
ing and secular coupling. In order to assess the strength of 
the possible interactions, we draw our initial eccentricity dis¬ 
tribution from the eccentricities observed in other planetary 
systems. Several proposals have been made for describing 
the exoplanet eccentricity distribution (e.g. Shen & Turner 
2008) but including planets at large orbital periods and sin¬ 
gle planets may overestimate the eccentricities for planets 
in compact multiple systems. Wang & Ford (2011) estimate 
the eccentricities for single planets with short orbital peri¬ 
ods but do not include multiples. Limbach & Turner (2014), 
on the other hand, derive the eccentricity distribution as a 
function of multiplicity, but include systems on large scales. 
We opt to draw our eccentricities from the ‘short period’ 
beta distribution ht by Kipping (2013) to the observed ex¬ 
oplanet data with orbital periods < 382 days. This draws 
values from the cumulative probability 


P(e) 


r(3.967) -0.303 _ 

r(0.697)F(3.27) '' 


( 10 ) 


We adopt this as representative of the potential initial con¬ 
ditions in the 55 Cancri system as it includes the bulk of the 
55 Cancri planets (except for 55 Cane d) within its period 
(< 387 days) threshold. The systems are then integrated for 
1 Myr. Our implementation of tidal evolution in Mercury6 
uses the force law used by Hut (1981), to damp eccentric¬ 
ity. As discussed in appendix we calibrate the dissipation 
in order to produce an eccentricity damping timescale of 
1.3 X 10® years for an 8M^ planet at a=0.033 AU. For a 
planetary radius of R — 2R^, this is equivalent to Q' = 13. 

Of the 100 planetary systems integrated, 33 underwent 
dynamical instability, usually after less than 0.1 Myr (where 
we define this as undergoing planetary scattering and colli¬ 
sions which removed at least one and often several planets 
from their observed orbits - be it by ejection or collision with 
the star or another planet). This indicates that even the re¬ 
duced eccentricities of Kipping’s ‘short period’ distribution 
are sometimes too large for a system as closely packed as the 
55 Cancri system. Of the 67 systems whose planets survive 
to the present, 38 reach semi-major axes < 0.03 AU within 
1 Myr, while 29 do not. Those that pass within 0.03 AU 
usually lie well inside 0.03 AU, as shown in the right-hand 
panel of Figure [T] This is a consequence of the eccentricity 
jump due to the secular resonance and the subsequent cir- 
cularisation of the orbit. Figure [7| also shows the resulting 
inclination distribution of the 38 systems with final semi¬ 
major axis < 0.03AU. We find inclinations can be pumped 
up to almost 60°, in agreement with the simple estimates 
of earlier sections. It should be noted that the inclinations 
shown in Figure [7] are the amplitudes of the secular oscilla¬ 
tions, not the instantaneous values, which is why the initial 
conditions (dotted histogram) are larger than the seed dis¬ 
tribution (dashed histogram) - the secular forcing from the 
other planets drives the inclination of 55 Cane e up even if 
the initial value is small. 

Of the 29 stable systems that do not penetrate inside 
0.03 AU within 1 Myr, it is possible that some may just 
have insufficient amplitude in the relevant modes to sec¬ 
ularly pump the eccentricity of 55 Cane e to values that 
promote tidal evolution. However, the distribution of initial 
eccentricities for 55 Cane b, which is the principal driver of 
the secular pumping, is quite similar to those cases which 
do migrate, as shown in Figure [S] This suggests that many 
of the systems have sufficient amplitude, but have had their 
eigenvalues shifted enough by the 3:1 resonance that it lies 
exterior to the starting location of our simulation. 

To test this hypothesis, we repeat the integration of 10 
of the systems which do not undergo sufficient evolution, 
with all of the same parameters except that the initial semi¬ 
major axis of 55 Cane e is taken to be 0.045 AU instead. 
This is to place it closer to the region destabilized in the 
case shown in Figure [6l Of these, two exhibited little evolu¬ 
tion, but in another five, 55 Cane e experienced such strong 
eccentricity excitation that it collided with the star with lit¬ 
tle semi-major axis evolution or inclination excitation. In 
the remaining three, the inner planet achieved a 55 Cane e 
analogue orbit, with a circularised orbit with a < 0.03 AU, 
and maximum inclinations of 8.2, 12.7 and 18.2°. Thus, in 
many of these systems, the inner planet experiences rapid 
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Figure 7. The right-hand panel shows the distribution of semi¬ 
major axis for the 38 systems whose evolution puts them interior 
to 0.03 AU. The observed location of 55 Cane e is shown as a ver¬ 
tical dashed line. The left-hand panel shows the final inclinations 
of these same systems (solid histogram). The dashed histogram is 
the distribution from which the initial planetary osculating incli¬ 
nations were drawn. Because of secular forcing, the actual initial 
proper inclinations are larger, and are shown by the dotted his¬ 
togram. The solid point with error bars indicates the inclination 
of 55 Cane d measured by McArthur et al. (2004), relative to an 
edge-on orbit, such as that of 55 Cane e. 

eccentricity growth, to the point that it rises to a catas¬ 
trophic level before tides even play a role. This is likely due 
to the additional forcing discussed in appendix lUl due to the 
3:1 resonance between 55 Cane b and c. 

Figures [9] and show two examples of the kind of evo¬ 
lution that might produce a misaligned 55 Cancri e. In Fig¬ 
ure we show an evolution from 0.033 AU which follows 
the path predicted by classical secular theory, in which a 
steady inward migration is punctuated with sudden jumps 
in eccentricity and inclination substantially in excess of the 
normal secular oscillations, resulting in a final inclination of 
almost 60° with respect to the original orbital plane. Fig¬ 
ure [ 10 ] shows the evolution in the case where the eccentricity 
is pumped up to large values even at a starting location of 
0.045 AU, which then drives tidal migration and again ex¬ 
cites inclination as it passes through the unshifted inclina¬ 
tion resonance. 

3.3 Effect of the Companion 

The presence of an M-dwarf companion at > 1000 AU can, 
in principle, induce precession of the planetary system (In- 
nanen et al. 1997; Kaib et al. 2011; Boue’ & Fabrycky 2014) 
if the orbit is sufficiently eccentric and inclined. However, 
the precession timescales are long compared to the time it 
takes 55 Cane e to cross the resonance, and so the orbit 
tilting discussed here is insensitive to the influence of the 
external companion. We have verified this by repeating five 



Figure 8. The solid histogram shows the initial eccentricity of 
planet 55 Cane b in the case where the planet e evolves to semi¬ 
major axis < 0.03AU within 1 Myr. The dotted histogram shows 
the initial eccentricity in the case where the system remains stable 
but does not satisfy this criterion. The dashed histogram shows 
the initial conditions that lead to dynamical instability and shows 
that eccentricities > 0.25 for planet b usually lead to instability. 
The similarity of the dotted and solid histograms suggest that it 
is not simply a question of secular amplitudes that determines the 
tidal evolution of the system. These distributions are consistent 
with a Gaussian distribution of eccentricities with dispersion = 
0 . 12 . 

of the simulations which crossed the resonance but includ¬ 
ing a O. 26 M 0 companion with semi-major axis 1250 AU, 
eccentricity= 0.93 and inclination 115° as discussed by Kaib 
et al. (2011). In all cases the eccentricity and inclination 
excitation was reproduced, with rapid eccentricity pumping 
followed by circulation, and the final orbital tilt of 55 Cane e 
relative to the other planets. 


4 RESULTS 

Our principal goal in this paper is to examine the effect of 
secular interactions coupled with tidal evolution on the con¬ 
figuration of the 55 Cancri planetary system. The generic 
behaviour of multi-planet systems under the influence of 
tidal evolution is that the eccentricities of several planets 
can be reduced, while the bulk of the semi-major axis evolu¬ 
tion is borne by the innermost planet (Wu & Goldreich 2002; 
Greenberg & van Laerhoven 2011; Hansen & Murray 2014). 
Such behaviour in this system would imply that 55 Cane e 
began with a larger semi-major axis and migrated to the one 
observed today. 

If 55 Cane e began with a semi-major axis > 0.033AU, 
it will have crossed multiple secular resonances, dominated 
by the pair of planets 55 Cane b and 55 Cane c, which could 
have excited both substantial eccentricities and inclinations. 
We have demonstrated these effects using both classical sec- 
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Figure 9. The upper panel shows the evolution of the semi-major 
axis, the middle panel the evolution of the eccentricity, and the 
lower panel the evolution of the inclination. The case shown here 
proceeds along the normal evolutionary path, with both eccen¬ 
tricity and inclination resonances occurring near 0.03 AU. 

ular theory, analytic estimates and direct numerical integra¬ 
tions. These suggest that inclinations of 30 — 40° between 
55 Cane e and the other planets in the system are quite pos¬ 
sible for initial conditions that draw planetary eccentricities 
from a Gaussian with dispersion (Je = 0.12 (an approximate 
fit to the dynamically stable portion of our simulations - 
see Figure [Hll and mutual inclinations from a Gaussian dis¬ 
tribution with a (Ti = 2° dispersion (consistent with the 
flat planetary systems observed by the Kepler satellite). The 
substantial amplification is a consequence of the large mass 
ratio between 55 Cane e and 55 Cane b, which causes any 
coupling between the planets to be reflected disproportion¬ 
ately in the smaller of the pair. The median eccentricity of 
the input planets, once the dynamically unstable systems are 
removed, is ~ 0.07, which is similar to the value expected 
for a five planet system (~ 0.085) from the power law fit of 
Limbach & Turner (2014). 

The ability to generate a substantial inclination offset 
from the orbital plane of the original planetary system of¬ 
fers an explanation for the tension between the astrometric 
measurement of 55 Cane d’s inclination by McArthur et al. 
(2004) and the transitting nature of 55 Cane e (Winn et 
al. 2011; Demory et al. 2011). One observation potentially 
in conflict with this scenario is the statistically marginal 
detection of an extended H atmosphere of 55 Cane b by 
Ehrenreich et al. (2012). If confirmed, this would suggest a 
different configuration for the system, with an approximate 
alignment of 55 Cane e and b. 

4.1 Planetary Heating 

The passage through the eccentricity resonance causes a 
rapid increase in the eccentricity of 55 Cane e and a cor¬ 



T (yrs) 

Figure 10. The upper panel shows the evolution of the semi¬ 
major axis, the middle panel the evolution of the eccentricity, 
and the lower panel the evolution of the inclination. In this case, 
the eccentricity resonance is shifted outwards, to ~ 0.043 AU. 

responding rapid decrease in the periastron distance. The 
strong radial dependance of the tidal force means that this 
leads to a sharp increase in the rate of tidal dissipation 
and tidal heating, which in turn may substantially influence 
the structure of the planet. In particular, we note that the 
amount of energy available could potentially have removed 
a substantial amount of Hydrogen from 55 Cane e. 

To estimate the impact of the heating on the structure 
of the putative planet, we note that the final planetary semi¬ 
major axis is roughly half that of the secular resonance, so 
that the binding energy dissipated as a result of the res¬ 
onance crossing is > 10“*^ ergs. This is comparable to the 
gravitational binding energy of a Neptune mass planet, and 
so can substantially influence the structure of the planet. 

The 55 Cancri system already contains several planets, 
such as 55 Cane c and 55 Cane f, which have masses simi¬ 
lar to Neptune or Uranus. Thus, 55 Cane e may potentially 
represent the rocky core of a thermally altered version of 
this class of planet. To ascertain the likelihood of this, we 
consider initial parameters for 55 Cane e, of a mass of 2 OM 0 
and a radius of 4i?0. The rate at which this energy is dissi¬ 
pated will depend on the rate of tidal dissipation, which we 
take to be characterised by Q' ~ 10"^, as in the case of Nep¬ 
tune (Banfield & Murray 1992; Zhang & Hamilton 2008) As¬ 
suming an average eccentricity ~ 0.5 during this period, and 
the constant time-lag model (Hut 1981) for tidal dissipation, 
we find that the binding energy is dissipated over a period 
~ 10® years, yielding an average luminosity ~ 10 -^Lq. 

The effective temperature of a 4i?0 planet, radiating at 
is ~ 3000K, comparable to the hottest estimated 
equilibrium temperatures of irradiated giant planets (e.g. 
Konacki et al. 2003; Hebb et al. 2009; Gillon et al. 2014), and 
~ 40% larger than the equilibrium temperature maintained 
by irradiation for this planet. Nevertheless, as in the case 
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Figure 11. The upper panel shows the evolution of both the 
semi-major axis (upper, smooth, curve) and the periastron (lower 
curve) for a 20Mq planet driven to high eccentricity by the cross¬ 
ing of a secular resonance. The horizontal dashed line indicates 
the semi-major axis at which a 20Mq, 4R0 planet will fill its 
Roche lobe. The lower panel shows the corresponding eccentric¬ 
ity excitation. 


of hot Jupiters on these scales, mass loss by Jeans escape 
is still not sufficient to substantially alter the planet mass, 
because A ~ GMuip/RkT ~ 12 for such a planet with this 
effective temperature. 

However, dissipation does not only increase the temper¬ 
ature of the planet. We have previously reviewed the impact 
of strong tidal heating on gaseous giant planets (Hansen 
2012) and found that the principal repository of the energy 
at high levels of dissipation is in the binding energy of the 
planet (i.e. at the level of dissipation relevant here, the ra¬ 
dius of the planet increases, and central temperature can 
even drop). For the parameters of relevance here, the radius 
need only expand to ~ 6.5i?0 to begin Roche lobe overflow 
at periastron and initiate mass loss. 

To demonstrate this, we have repeated one of our sim¬ 
ulations but replaced 55 Cane e with a 2 OM 0 object, and 
adjusted the level of dissipation to mimic a value of Q' ~ 10^ 
as appropriate for Neptune. Assuming that the radius of the 
planet is now 4f?0, this requires weakening the numerical 
dissipation by only a factor of 60, because Q' scales with 
the fifth power of the planetary radius. The resulting evo¬ 
lution is shown in Figure 1111 We see that the character of 
the evolution is the same, with excitation to high value of 
eccentricity and a sharp drop in semi-major axis, although 
a little less steep, as befits the weaker dissipation. We also 
show the evolution of the periastron during this period and 
note that it skirts the edge of Roche-lobe overflow even for 
an un-infiated radius of 4i?0. Thus, even a small amount of 
radius inflation due to the tidal dissipation will drive this 
planet into a mass loss event. 

Thus, it is possible that the present mass and radius 


of 55 Cane e are the result of mass loss, probably through 
Roche lobe overflow, induced by the strong tidal heating 
generated by the passage through the secular eccentricity 
resonance. If this is the case, 55 Cane e may originally have 
been a Neptune-mass planet, more similar to other known 
members of the system, such as 55 Cane c and 55 Cane f, 
suggesting migration as a chain of giant planets and failed 
cores. Residual heating may still make a contribution to the 
amplitude of the secondary eclipse (Bolmont et al. 2013) if 
the eccentricity is in the range ~ 0.001 — 0.01. 

4.2 Other Systems 

The results above suggest that the 55 Cancri system may 
exhibit a substantial inclination dispersion, but this is also 
a relatively unusual planetary system, as most short period 
Jupiters do not reside in systems of non-hierachical high 
multiplicity (Wright et al. 2009; Steffen et al. 2012). Thus, 
it is of interest to enquire as to the frequency of such secu¬ 
lar resonances in compact planetary systems. Table 2 shows 
examples of several systems which may possess such a reso¬ 
nance. 

For systems detected via radial velocities, we examine 
those systems which possess three or more planets inside 
0.5 AU, at least two of which lie outside 0.05 AU. In ad¬ 
dition to 55 Cancri, this yields 61 Vir (Vogt et al. 2010), 
HD 40307 (Mayor et al. 2009) and HD 10180 (Lovis et al. 
2011) around solar-type stars, and GJ 876 (Rivera et al. 
2010) around an M star. In each case, we consider the secu¬ 
lar architecture of the system given the nominal semi-major 
axes and masses, and examine how this changes as we vary 
the semi-major axis of the innermost planet. This accounts 
for possible motion due to tidal evolution, which is reflected 
primarily in the evolution of the innermost planet for sec¬ 
ularly coupled systems (Hansen & Murray 2014). However, 
some care is necessary as several of these systems contain 
uncertain detections. For HD 40307, we have adopted the 
more conservative configuration because the more uncertain 
planets (Tuomi et al. 2013) reside further out and including 
them turns out to not significantly change the results. In 
the case of HD 10180, we include the uncertain innermost 
planet as a potential analogue to 55 Cane e. We note also 
that more strongly resonant systems, like GJ 876, may ex¬ 
hibit more complicated behaviour due to untreated resonant 
effects. 

With these caveats, we find that each of the examined 
systems possess at least one secular resonance in both incli¬ 
nation and eccentricity in the region through which a planet 
would pass due to tidal evolution. In some cases, such as the 
61 Vir and HD 10180 systems, the inner planet would ap¬ 
pear to actually be near the location of a secular resonance 
at present, offering a potential explanation for non-zero ec¬ 
centricities despite proximity to the star yielding strong tidal 
effects. 

In the case of systems detected via transits, passage 
through an inclination resonance is likely to prevent a planet 
from being observed in transit, assuming the natal planetary 
system is edge-on in the first place. Such excitations may 
explain inner holes in otherwise densely packed systems ob¬ 
served by Kepler. An example is the Kepler-11 system (Lis- 
sauer et al. 2011b), which shows 6 planets, five of which 
are between 0.09-0.3 AU, but none interior to 0.09 AU. 
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Table 2. Orbital periods of possible secular resonances in exoplanetary systems, where Pe indicates ec¬ 
centricity resonances and indicates inclination resonances. Numbers in parentheses are caused by first 
order mean motion resonances. The final column for the RV planets shows the observed orbital period of 
the currently known exoplanet. 


System 

Pe 

(days) 

Pi 

(days) 

P 

(days) 

RV Systems, assuming inner planet migration 

p 55 Cancri 

1.66 

2.0 

0.74 

GJ 876 

7.1, (14.7) 

7.3, 9.8, 17.0 

1.9 

61 Vir 

4.2, (18.5) 

6.8 

4.2 

HD 40307b 

0.9, (4.7) 

to 

to 

CO 

4.2 

HD 10180 

1.2, (3.0) 

0.8, 1.6 

1.2 

Transitting Systems, assuming an undetected IM 0 inner planet 

KOI -94 

0.9, (1.87) & 2.1 

1.2 & 2.1 


Kepler-11 

2.4, (5.1), 8.0 

1.3, 6.1, 8.0 


Kepler-30 

(13.9), 14.8 

14.8 


Kepler-9 

1.0, 2.5 

7.5 


Kepler-51 






Figure 12. The curves show the change in the secular inclination 
eigenfrequencies of a hypothetical planetary system constructed 
by adding an interior IM^ planet to the known planetary sys¬ 
tem Kepler-11. As the semi-major axis of the additional planet is 
varied, we see that there is an avoided crossing at 0.065 AU. 

If we examine the secular architecture of this system as¬ 
suming the presence of a seventh, IM^ planet located in 
the interior holfEl, we find avoided crossings in both eccen¬ 
tricity and inclination near semi-major axes ~ 0.065 AU, 
and other secular inclination resonances at 0.0233 AU and 
0.0763 AU. Thus, it is possible that Kepler-11 contains a 

^ Such a planet would yield transit timing variations of amplitude 
< 1 minute for Kepler-llb (Agol et al. 2005), well below the 
current timing accuracy (Lissauer et al. 2013), as long as it lies 
interior to the 2:1 mean motion resonance. 


low mass planet in the observed inner ’hole’, but which is 
inclined relative to the outer planets. The inclination could 
be substantial if the current period is < 8 days. The Kepler 
planet sample contains a large number of multiple planet 
systems with potential for secular resonances. KOI-94 (Weiss 
et al. 2013) is interesting because of the presence of a Sat¬ 
urn mass planet in a similar orbit to 55 Cancri. The inner¬ 
most planet lies just outside a pair of secular resonances at 
3.0 days (eccentricity) and 3.4 days (inclination) and so is 
consistent with a lack of inclination excitation. If we allow 
for the existence of a IM® planet interior to this, eccentric¬ 
ity resonances are found at 0.9 and 2.1 days (and a mean 
motion resonance at 1.9 days) and inclination resonances at 
1.2 and 2.1 days. Thus, a planet that has tidally migrated 
to orbital periods <2 days in this system could have had its 
inclination pumped up and have avoided transit. Another 
compact multiplanet system containing a massive planet is 
Kepler-30 (Fabrycky et al. 2012). However, the only reso¬ 
nances encountered by a tidally migrating IM^ inner planet 
lie at 14.8 days (both inclination and eccentricity) and at 
13.9 days (a mean motion resonance). These may be too 
far out to be strongly affected by tides. Kepler-9 (Holman 
et al. 2010), shows two Saturn-class planets and an inte¬ 
rior low mass planet, that we will assume to be IM® again. 
Since this inner planet does indeed transit, we must assume 
little inclination excitation has taken place in this system. 
Indeed, the sole inclination resonance for this configuration 
lies at 7.5 days, so that the planet would not experience any 
excitation if it started the migration interior to this point. 
Eccentricity resonances are found at 1.0 and 2.5 days, brack¬ 
eting the observed period. Kepler-51 (Steffen et al. 2013) is 
an example of a compact Kepler planetary system with a 
large inner hole but no Jovian-class planets. In this case, an 
inner IM® planet experiences no secular resonances. There 
is a frequency commensurability at 2.9 days, but no avoided 
crossing because the modes are almost completely decoupled 
(one is driven primarily by the relativistic precession). 

These results are necessarily incomplete because our 
knowledge of the planetary systems could also be incom¬ 
plete, but they illustrate the potential impact of secular res- 
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onances on the internal planetary configurations of multi¬ 
planet systems. They may be less important for lower mass 
systems where the precessions are weaker and relativistic 
effects can dominate at larger distances, but could gener¬ 
ate substantial inclinations in the cases where Jovian-class 
planets are present on small scales. Alternatively, the effect 
of stellar spin-down from an initially rapid rotation may, in 
some cases, also induce secular resonant crossings, as was 
once proposed for the eccentricity and inclination of Mer¬ 
cury (Ward, Colombo & Franklin 1976). 


5 CONCLUSIONS 

We find that it is possible to explain the conflicting obser¬ 
vations that 55 Cane e transits, while 55 Cane d is inclined 
by ~ 53°, in terms of tidal decay of the 55 Cane e orbit 
that causes it to cross multiple secular resonances, which 
drive the planet to large eccentricities and inclinations, van 
Laerhoven & Greenberg (2012) have previously cast doubt 
on whether the system can have undergone tidal evolution 
because not all of the secular modes coupled to the inner 
planet are completely damped. However, a system need not 
have reached it’s final state of a ‘fixed point’ configuration, if 
the evolutionary timescale is not too short. Indeed, Hansen 

6 Murray (2014) found many model systems which were 
only partially damped, and the large amplitude librations 
about secular alignment found by Nelson et al. are char¬ 
acteristic of this stage of evolution. Furthermore, it worth 
noting that the updated parameters of Nelson et al. do dif¬ 
fer in important ways from the solutions analysed by van 
Laerhoven & Greenberg (2012). The pericenter alignment 
of the outer, decoupled, pair of planets, which appeared to 
be coincidental, is no longer supported by the data. The re¬ 
maining eccentricity of 55 Cane c suggests that this system 
is indeed still capable of further tidal evolution, although the 
presence of a tentative signal (Nelson et al. 2014) between 
planets 55 Cane c and 55 Cane f may hint at additional 
influences that can masquerade as an eccentricity. 

The evolution through secular resonance increases the 
eccentricity as well as the inclination for this planet, and can 
result in a substantial amount of heating within 55 Cane e. 
This has potentially important consequences for the use of 
55 Cane e as one of the highest signal-to-noise examples of a 
super-Earth planet. We have shown that it is possible that 
55 Cane e could have lost mass through Roche lobe over¬ 
flow if it expanded due to tidal heating, and so 55 Cane e 
may represent a class of disrupted planetary cores that 
have been occasionally invoked by several authors (Gu, Lin 
& Bodenheimer 2003; Jackson, Barnes & Greenberg 2009; 
Hansen 2012; Valsecchi, Rasio & Steffen 2014), rather than a 
rocky planet that assembled through planetesimal collisions. 
A more massive and volatile-rich initial mass would make 
55 Cane e comparable with the other sub-jovian planets in 
the system, and suggests a scenario in which 55 Cane b mi¬ 
grated inwards with a retinue of Neptune-mass attendants, 
one of which was subsequently heated and stripped to the 
planetary core we observe day. 

We have also shown that many compact multiple 
planet systems exhibit possible secular resonances interior 
to 0.1 AU, suggesting that the phenomenon of moderate-to- 
highly inclined inner planets may be quite common. Such 


inclination pumping could potentially explain the observed 
lower deficit of multiple transitting planets at very short pe¬ 
riods in the Kepler sample (Steffen & Farr 2013). Similarly, 
pumping of eccentricity and subsequent tidal circulation is 
a potential mechanism for moving planets to very short or¬ 
bital periods and promoting their disintegration (Rappaport 
et al. 2012, 2014; Jackson et al. 2013). 

The analysis presented here is predicated on the HST 
astrometric detection by McArthur et al. (2004), which 
places some constraints but is based on a limited data span. 
Further astrometric observations of this system, either with 
the HST Fine Guidance sensor or the GAIA spacecraft, 
would provide a welcome improvement on this constraint, 
and could help to refine the parameters. A substantial mu¬ 
tual inclination between 55 Cane e and the rest of the planets 
could also be tested with transit monitoring, to search for 
precession of the orbit. 

This research has made use of NASA’s Astrophysics 
Data System and of the the NASA Exoplanet Archive, which 
is operated by the California Institute of Technology, under 
contract with the National Aeronautics and Space Admin¬ 
istration under the Exoplanet Exploration Program. The 
authors also acknowledge the contributions from a useful 
referee report. 
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APPENDIX A: ECCENTRICITY AND 
INCLINATION EXCITATION DUE TO 
SECULAR RESONANCE CROSSING 


We consider an extended version of the Hamiltonian dis¬ 
cussed in Minton & Malhotra (2011), 


H 


Re + Hi (Al) 

—goJ — hoZ + eVn cos(a;p — cu) -b &\f2Z cos(flp — fl) 


where J = \/a(l — v/l ~ e^)) Z = o(l — e^)(l — cost), 
go and ho are the precession frequencies for the longitude 
of periastron tj, and longitude of the ascending node Q re¬ 
spectively. The parameters e and S represent the strength 
of the coupling to the secular modes, whose orientations are 
represented by uip, and flp. 

We express these quantities as 


go 

ho 


2,(1) / N 


SGrric 


3 

m 


3 

= J OGnica 4 %(aj )E, 


C n "ij 2,(1)/ w 


(A2) 

(A3) 

(A4) 

(AS) 


where Ej and Ij are the eigenvector contributions from each 
planet to the mode under discussion, and go includes the 
effect of relativistic precession to lowest order. 

The approach of MMll is to perform a pair of canonical 
transformations to variables x and y, represented by a; = 
'JTj cos (6 and y = 'JTj sin((), where <t> = CJp{t) — oj. In our 
case, it is oj, rather than ojp, that is a function of time, but 
the formalism remains the same. However, we must now 
incorporate changes in e concommitant with the changes in 
go, since both are affected by the radial migration of the 
planet. 

Following the formalism of MMll, we find that the 
Hamiltonian equations of motion yield 


X = goyt (A6) 

y = -e - {gox + k)t (A7) 

which is the same as equations (9) and (10) of MMll, except 
for the £ term. 

The method of solution is similar, whereby one can find 
a homogeneous solution in terms of trigonometric functions 
of goT, and then substitute back to find a solution to the 
inhomogeneous version by solving for the multiplicative con¬ 
stants as functions of time. Thus, we assume 

x{t) = A(t) cos goT -b B{t) sin goT (A8) 

y{t) = —A{t) sin goT + B{t) cos goT (A9) 


MMll solved for A and B in terms of the Fresnel inte¬ 
grals 


S{t) = ( sin z^dz 

Jo 

C{t) = [ cos z^dz. 

Jo 


(AlO) 

(All) 
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The presence of e term introduces a slight modification to 
the solution, yielding 


A{t) = + 

l5o| 


B{t) = - 






(A12) 

(A13) 


where t' = ^J\go\t, and S' and C are modified versions of 
the Fresnel integrals, such that 


and 


s'it) = f z sin z'^dz 
Jo 

H) ~ [ z cos z^dz. 
Jo 


For our purposes, we care primarily about the prop¬ 
erties of the solution as t' goes from —oo to -|-oo (i.e. the 
transformation of the orbit as it passes well beyond the res¬ 
onance). Thus, we can make use of the identity 


n \ n J 


^i7r(l + m)/2n 


to evaluate S'(oo) = C(oo) — S{—t) = —S{t), 

C{—t) = —C(t) (as in MMll), as well as S'{oo) = 1/2, 
C’'(oo) = 0, and S'{-t) = S'{t), C'{-t) = C'{t). So, we 
find that the evenness of the functions S' and C' mean that 
the contributions of the e term to the asymptotic behaviour 
average out and we recover the same result as MMll. 

Thus, our final expressions are the same as those from 
MMll, 


Jf — Ji ■ 


+ e 


27r Ji 


2|5o| V ISol 

with the equivalent expresson for the inclination 

„ ry „ / 27rZi 

Zf — Zi -\- —p—j- -p—p cos fl. 

2\ho\ Y \ho\ 

The expressions for e and S are as above, and the ex¬ 
pressions for go and ho are 


1 V—' nr, 2 

3 


dgo 1 da 

dt a dt 

and 

dho 1 da 1 

—— = — -n- > — -ctl'Ji 

dt a dt 8 rric ^ 


15 Grric 


c^a 




(A14) 


(A15) 


where jj = -I- Saj 

For the 55 Cancri system, the relevant frequencies and 
couplings are all dominated by 55 Cane b, with contributions 
from 55 Cane c and relativity contributing at the one percent 
level. The eigenfunction for the mode that drives the excita¬ 
tion through the resonance actually contains a stronger con¬ 
tribution (see Figure [1} from 55 Cane c (by a factor ~ 5), 
but the greater mass and proximity of 55 Cane b means 
that the functions e and 5 still only contain a ~ 5% con¬ 
tribution from 55 Cane c. Thus, in our discussions, we will 
refer primarily to excitation by interaction with 55 Cane b. 




0.006 


0.004 


0.002 


0 



0 106 2x106 3x106 4x106 

T (days) 


Figure Bl. The points show the accumulated deviation between 
the numerically integrated argument of periastron and the value 
calculated analytically, for a planet of semi-major axis 0.05 AU, 
and eccentricity e = 0.2056. Only a single planet was included 
and no tidal damping was included. 


APPENDIX B: CALIBRATION OF THE 
RELATIVISTIC PRECESSION AND 
NUMERICAL TIDAL DISSIPATION 

The Mercury code (Chambers 1999) does not include a pre¬ 
scription for relativistic precession or tidal damping, so we 
incorporate an extra user-defined force as follows. The rel¬ 
ativistic precession term is standard but can sometimes re¬ 
quire special handling for proper implementation (e.g. Saha 
& Tremaine 1994). However, our interests are on sufficiently 
small scales that our timesteps are chosen small enough to 
properly resolve the precession anyway, and a direct imple¬ 
mentation with the MVS symplectic algorithm yields an¬ 
swers of sufficient accuracy. Figure IBII shows the accumu¬ 
lation of the difference between the numerically calculated 
relativistic precession using Mercury and the expected an¬ 
alytic value (calculated for a single planet at 0.05 AU and 
with no tidal effects included). The drift of the location of 
periastron is ~ 5 x 10~® °.yr~^, well below any of the pre¬ 
cession frequencies of interest in this system (see Figure [T]). 

We have previously (Hansen 2010) calibrated tidal in¬ 
teractions based on the constant time-lag model of Hut 
(1981) and Eggleton, Kiseleva & Hut (1998), and so we 
adopt the force law from Hut (1981). The component that 
damps eccentricity is that which acts to dissipate the radial 
component of the planetary velocity, namely 



where r and r are the heliocentric radius and radial veloc¬ 
ity vectors. In Hut’s formula, the constant a — dCMH^Rp, 
where k is the apsidal constant, r a timescale that encodes 
the strength of dissipation, and Rp is the radius of the planet 
(in which the tides are being dissipated). For the purposes 
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Figure B2. The solid line shows the evolution of an isolated 
55 Cane e analogue, under the influence of our tidal dissipation 
prescription. The plot actually contains a dashed line representing 
e = 0.5exp(—4/1.3 X lO^yrs), but it is hard to see because the 
fit is excellent. The dotted line relates the eccentricity to the 
semi-major axis, given on the upper axis, and demonstrates the 
conservation of angular momentum during tidal decay. 


of calibrating the simulations a is a single number, whose 
physical relevance will depend on the relationship between 
the various physical quantities of which it is composed. In 
the low eccentricity limit, the resulting evolution of the ec¬ 
centricity is an exponential damping, whose characteristic 
timescale is determined by a. 

Figure IB2I shows the evolution for a 55 Cane e ana¬ 
logue, begun at 0.033 AU, with eccentricity=0.5 with the 
same calibration for a as the simulations described in the 
main text. This is very well fit by an exponential decay on 
a timescale of 1.3 x 10® years. In the formalism of Hansen 
(2010), and assuming a planet of mass 8M® and radius 27?®, 
this corresponds to ap ~ 2.5 x g~^cm~^s~^. This is a 
considerably stronger level of bulk dissipation than found 
for giant planets in that paper, but such a difference is to 
be expected for a terrestrial-class planet, and indeed, it is 
only a factor of 10 larger than the nominal value used by 
Bolmont et al (2013), based on an extrapolation from Earth 
dissipation levels. If cast in terms of the traditional tidal Q, 
it corresponds to Qp = 13, using the expression from Jack- 
son, Greenberg & Barnes (2008), which is also comparable 
with an earth-like dissipation. 

If we adopt a Neptune mass and radius of 20M® and 
47?®, this would correspond to Qp ~ 150, which is a little low 
compared to that observed for Uranus & Neptune (Banfield 
& Murray 1992; Zhang & Hamilton 2008). 

Figure IB2I also shows the trajectory of the planet evo¬ 
lution in semi-major axis as a function of eccentricity. The 
curve traces out that expected due to the conservation of an¬ 
gular momentum during isolated tidal evolution. It is also 
worth noting that the eccentricity properly asymptotes to 


zero at late times. This can be compared with the late-time 
behaviour shown in the full runs in Figure [9l and llOl which 
show some fluctations. The fact that they do not occur here 
demonstrates that the effect is a real one, due to residual 
secular perturbations from the other planets, and not due 
to numerical inaccuracies. 


APPENDIX C: RESONANT INTERFERENCE 
WITH THE SECULAR INTERACTIONS 

The numerical experiments demonstrate that the planets 55 
Cane b and c are close enough to the 3:1 resonance that 
their resonant interactions may sometimes distort the secu¬ 
lar structure of their interactions with 55 Cane e. Figures [5] 
and [6] showed how this can manifest itself in the Lagrangian 
forcing of test particles. A more direct demonstration can be 
found in the Fourier decomposition of an integration of the 
full planetary system. We have performed a fourier trans¬ 
form of the time series of the eccentricity of 55 Cane b for 
each of the two cases shown in those figures, 

N-l 

Ck = J2 (Cl) 

1=0 

and calculated the power spectral density 

P{f) = ^[\Ck\^ + \CN-k\^], (C2) 

for a time series with a 10 day timestep and N = 2^^. 

Figure El] shows the power spectrum P{f) for the case 
shown in Figure [S] which conforms to the standard sec¬ 
ular evolution we expect. The power spectrum shows the 
strongest power at the expected secular frequency gi (the 
time series was not long enough to capture the longer sec¬ 
ular period) with additional power in the neighbourhood 
of the circulation frequency of the near resonance. There is 
also a minimum in the power between these two values, sug¬ 
gesting a clear separation between the resonant and secular 
phenomena. 

For a different set of initial eccentricities and longitudes 
of periastron, corresponding to Figure [b] we find the power 
spectrum shown in Figure [C2| In this case the secular fre¬ 
quency is shifted to higher values relative to the classical 
value gi, and more power is found between this frequency 
and the higher circulation frequency, indicating a stronger 
interaction. The shift in frequency of the secular mode is 
consistent with the change in the location of the secular res¬ 
onance observed in Figure [5| 

Could this simply be the case of a set of initial condi¬ 
tions that put almost no power into the mode with frequency 
gi, rather than a shift driven by proximity to resonance? Fig¬ 
ure [C^ shows the power spectrum of a simulation run with 
exactly the same parameters as in the case of Figure IC21 
except that the semi-major axis of 55 Cane c was moved 
from 0.237 AU to 0.24 AU. We see that this change in 1% is 
sufficient to restore the mode frequency to it’s expected clas¬ 
sical value, indicating that the observed change in frequency 
is directly linked to the proximity to resonance. 

To gain a qualitative understanding of this phe¬ 
nomenon, let us collect the relevant secular and resonant 
terms, up to second order in eccentricity, for the pair 
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Table Cl. Constants that define the averaged secular and resonant interactions of the 55 Cane b and 


55 Cane c planet pair. 


Cl 

0.263 

Ci 

2.996 

C 2 

0.139 

Cs 

-6.133 

C 3 

-2.149 

Ce 

1.197 



100 1000 104 10 ^ 10 ® 


1/f (days) 


Figure Cl. P{f) is the power spectral density of the eccentric¬ 
ity of 55 Cane b, plotted against the corresponding period in the 
time series. The vertical dotted line shows the period associated 
with the expected circulation frequency 3nc — rih of the planet 
pair 55 Cane b and 55 Cane c. The vertical dashed line is the 
expected secular frequency gi based on classical secular pertur¬ 
bation theory. 


55 Cane b and 55 Cane e. The disturbing funetions, av¬ 
eraged over short period terms, are 

< Rb >= 8.73 X 10"'‘ [Co + Cl {el + el) + 

Csetec cos {oje — Ojb) + CacI cos (0 — 2ujb) 

+Csebec cos {O — LUb — uic) + Ceel cos (0 — 2ajc)] (C3) 
< Rc 4.133 X 10 ^ [Co “t“ Cl {eb 4“ e^) ~\~ 

CsebCo cos {uic — Ojb) + C 4 e 6 cos (0 — 20jb) 

-\-Cbebec cos {Q — Ojb — oJe) + Cjel cos (0 — 2cDc)] , (C4) 

where the Ci are given in Table ED and the difference be¬ 
tween Co and C 7 is due to the differing contributions to the 
3:1 resonance from the internal and external components 
of the indirect disturbing function. Far from resonance, the 
terms in the angle 0 = 3Ac — \b will average to zero and 
the system will exhibit the traditional secular oscillations. 
However, if the secular oscillations bring the planetary ec¬ 
centricities to low enough values, the resulting increase in 
the precession rate can sweep the system through resonance, 
inducing a brief libration. 

We integrate the Lagrange equations based on equa¬ 
tions (ESI) and (E4l) to describe the evolution of the 



l/f (days) 


Figure C2. P{f) is once again the power spectral density of 
the eccentricity of 55 Cane b, but now for a case in which the 
behaviour deviates significantly from the expectations of classi¬ 
cal secular theory. In particular, we note that the period of the 
dominant secular mode is now substantially shifted from the ver¬ 
tical dashed line. A comparison with the structure observed in 
Figure IC2I demonstrates that the power in the frequency range 
between the secular frequency and the circulation frequency is 
higher, indicating a stronger degree of interaction. 


55 Cane b/55 Cane c system, for a set of initial conditions 
that showed the resonant interaction during the full numer¬ 
ical integration. 

As Figure lCil demonstrates. the evolution that arises ex¬ 
hibits the characteristic oscillatory behaviour of secular evo¬ 
lution, with the planetary eccentricities undergoing quasi- 
sinusoidal oscillations with a range of frequencies. The most 
rapid oscillations are associated with the circulation of the 
resonant argument. However, we also see the occasional 
transient librations which act to change the phase of the 
oscillations and equivalently the amplitudes of the oscilla¬ 
tions. The inset in Figure ES] shows the effect of one of these 
librations, which acts to essentially apply a random phase 
kick to the overall oscillations. Figure ES] demonstrates the 
nature of the deviations from simple circulation during this 
interaction. The four panels show a sequence of trajectories 
in the space Xb = eb cos (0 — 2CJb) and yb = eb cos (0 — 2u}b). 
When the eccentricity gets small, the trajectory exhibits a 
reversal, which represents a transition to libration. A trajec¬ 
tory for a more lengthy integration is shown in Figure IC61 
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Figure C3. P{f) is once again the power spectral density of the 
eccentricity of 55 Cane b, but now for a case in which the semi¬ 
major axis of 55 Cane c has been moved out by 1%, to reduce the 
effects of resonant interactions. All other parameters, including 
initial eccentricities and phases, are the same as in Figure IC2I 
We see here that the distinction between the secular interactions 
and the circulation of the resonant angle is much clearer, and 
that the values correspond to the expectatios of classical secular 
theory. 


which demonstrates the presence of a critical eccentricity 
below which the circulation transitions to libration. 

To understand this, let us extract the resonant Hamil¬ 
tonian from the relevant part of < Ri, >, to study the dy¬ 
namics of the el resonance at low Cb- Normalising by the 
mass, and defining (pb ~ O — and Auu = ujb — uic, we get 
the Hamiltonian 

H = \CAel -b CsebCc cos Aui + Ceel cos 2Atj] cos <j>b 

+ [CbebCc sin Auj -|- Ceel sin 2Aa;] sin (f>b- (C5) 

This Hamiltonian has several equilibrium points, but 
some are for the case when 66 and Cc are of similar magni¬ 
tude, which is not our focus here. In the case when 66 <C Cc 
and LJb = uJc, we find two equilibrium points at ®6 = 0 and 

\Vb\ = (-Cs ± [cl - 4C'4C'6]'''") = 0.2176.. (C6) 

This value of yb defines a characteristic value of et at 
which the circulation of the resonant argument switches over 
to libration. Figure ICtI shows the curves of constant Hb for 
the case Cc = 0.12 (corresponding to the value during the 
resonance passage shown in the inset of Figure [C4|, which 
yields the critical value of 66 = 0.026 (shown as a dashed 
circle). When the secular oscillations bring 66 below this 
value, the increase in secular precession sweeps the system 
into resonance, and the trajectory only returns to simple 
circulation when the oscillations bring 66 to larger values 
again. The behaviour is similar when the periapses are not 
aligned, but the locations of the equilibrium points rotate 
in Xb and yb- The dramatic behaviour seen in Figures IC5I 



T 


Figure C4. The solid points show the eccentricity of 55 Cane b, 
integrated using the Lagrange equations and disturbing func¬ 
tions El and El. while the open circles show 55 Cane c. The 
curves show longer timescale secular oscillations allied to shorter 
timescale resonant oscillations. The inset shows a zoom of the 
passage through a brief libration in the angle <f> = 3Ac — Xb~ 2^6. 
with the arrow indicating the time of the libration. This corre¬ 
sponds to a change in the phase of the secular oscillation. 


and IC6I are thus manifestations of the passage through the 
resonance as the rate of secular precession increases while 66 
decreases. 

Similar considerations describe the behaviour of the 
el resonance, when Cc ^ 66. Thus, although the nominal 
timescales for secular oscillations and mean motion libra- 
tions are separated by almost two orders of magnitude, the 
excursions to low eccentricities in the secular oscillations are 
sufficient to enhance the precession rates to sweep the sys¬ 
tem through the resonance on an intermittent basis. The 
resulting interference with the secular oscillations appears 
to shift the frequencies sufficiently to make a measureable 
difference in the resonant forcing of some of our model sys¬ 
tems. However, it appears as though this subset has less 
likelihood of generating observable 55 Cane e analogues, as 
most of the inner planets in these systems end up hitting 
the star (§ 13.211 . 
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Figure C5. Panels a-d indicate a sequence of trajectories for 
the quantities and ^ they pass through an episode of tran¬ 
sient libration. The particular interaction show here is also the 
one responsible for the evolution during the inset in Figure [C4l 
In particular, panel (a) covers the time 1100—1160, (b) from 1160— 
1230, (c) from 1230-1300 and (d) shows the evolution from here 
through the next interaction, i.e. from 1300—1600. The dotted 
circles in each panel have a radius of 0.025 and represent the crit¬ 
ical eccentricity below which one expects an unstable equilbrium 
point to appear in this resonance, as discussed in the text. 


Figure C7. The curves show lines of constant based on equa¬ 
tion llcsl l. for the case Cc = 0.12. The dashed line indicates a circle 
of e}) = 0.026. The two solid points indicate the equilibrium points 
on these scales. The sign of H changes upon passage through the 
critical at fixed x;,, which drives the kinds of reversals of tra¬ 
jectory seen in Figure ED because the path has to change from 
one side of the origin to the other. 



- 0.04 - 0.02 0 0.02 0.04 


Figure C6. The heavy dashed curve represents the critical ec¬ 
centricity of 0.025, and the smaller dots show the trajectory of 
the system evolution, which is characterised by a simple circula¬ 
tion exterior to this curve but by a more complicated and chaotic 
evolution interior to this value, as the system exhibits a resonant 
interaction. 
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